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J> ' Abstract 

OO '. 

£f*l , I address the current status of dynamical decoupling techniques in 

<^> terms of required control resources and feasibility. Based on recent 

advances in both improving the theoretical design and assessing the 
control performance for specific noise models, I argue that significant 
progress may still be possible on the road of implementing decoupling 
under realistic constraints. 



1 Introduction 
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^ Achieving decoherence control has become a critical goal for a variety of 

applications in contemporary physics and engineering. Broadly speak- 
ing, the challenge is to preserve quantum coherence during both the 
storage and the manipulation of quantum states in systems which are 
■ unavoidably exposed to the influence of their surrounding environment 

I - therefore are subject to irreversible, open-system dynamics A 

number of practical and conceptual motivations have contributed to 
make this challenge increasingly timely and important in recent years. 
On one hand, ensuring a sufficiently high degree of control over both 
environmental and operational errors is a chief requirement for physical 
realizations of quantum information, shared to a lesser or greater extent 
by the large variety of device technologies which are nowadays being 
considered |3J0]. In a broader context, maintaining coherent quantum 
behavior in the presence of realistic noise sources is essential for guaran- 
teeing proper performance by any nanoscale process or device intended 
to operate in a quantum-limited regime [1] . On the other hand, from a 
more conceptual standpoint, the development and characterization of 
control-theoretic notions and tools for generic open-system evolutions 
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is a fascinating area per se, which remains only partially explored to 
date 00. 

As a result, the field of dccoherence control has rapidly grown 
recently, and several quantum stabilization strategies have been de- 
veloped, both specifically tailored for quantum information protec- 
tion and processing, and beyond. These include passive stabilization 
schemes inspired to quantum reservoir engineering jS], decoherence- 
free-subspace [§] and noiseless-subsystem coding JU] (including topo- 
logical approaches as well as active stabilization techniques 
based on quantum feedback control quantum error-correcting 
codes ^] EL an d dynamical decoupling methods j^BJ ^T|- While a 
discussion of the principles underlying these different approaches is be- 
yond the present purposes, 1 will focus here on critically reconsidering 
the control resources involved in the standard dynamical decoupling 
setting, and argue that substantial progress might be possible by both 
further improving the control design and by better incorporating avail- 
able knowledge about the relevant environmental interaction. 

2 Dynamical decoupling framework 

Dynamical decoupling methods for open quantum systems derive their 
basic physical intuition from coherent averaging techniques in high- 
resolution nuclear magnetic resonance (NMR) spectroscopy E| > 
decoupling being thereby understood as the selective removal of un- 
wanted contributions to the nuclear spin Hamiltonian via the appli- 
cation of suitable pulse sequences. In the current control-theoretic 
terminology, decoupling is usually interpreted in a broader sense than 
the original one, see [II [T3 E01 El E2 E3 El E3 EBJ for recent repre- 
sentative literature. A definition which I will adopt within the present 
discussion is to generally understand a decoupling scheme as an open- 
loop dynamical control protocol which relies on the repeated applica- 
tion of pulsed or switched controls drawn from a finite set in order 
to attain a desired control objective. In what follows, I will focus on 
the prototypical application of decoupling, in which case the objec- 
tive is the averaging of the environmental couplings responsible for 
dccoherence. In general, however, alternative control scenarios may be 
envisaged, notably the synthesis of controlled evolutions according to 
predetermined symmetry criteria |171 12 1| , and Hamiltonian quantum 
simulations [57] . 

The relevant decoupling setting may be pictorially schematized as 
in figure 1. The joint evolution of the target system S in interaction 
with the environment E is described by a total drift Hamiltonian of 
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Figure 1: Block-diagram of an open-loop controlled quantum system, S, in inter- 
action with an environment. E. No direct control action on the latter is assumed. 



the form 

H = H S ®I E +I S ®H E + H SE , H SE =^2s a (g)E a , (1) 

a 

where Hg and H E account for the isolated dynamics of the system 
and the environment, respectively, and the interaction term Hg E is re- 
sponsible for introducing unwanted decohcrence and dissipation effects 
in the reduced dynamics of S alone. The basic idea is to adjoin a ju- 
diciously designed controller, described by a classical time-dependent 
field H c (t) acting only on the system, in such a way that the result- 
ing controlled dynamics is described by an effective Hamiltonian H e f / 
which no longer contains mixing terms between S and E that is, 

H a ^ H Q + H c (t) such that H ef f = H s ® I E + I s ® H E , (2) 

for an appropriate, possibly modified, system Hamiltonian H$. 

Decoupling protocols are most conveniently constructed by directly 
looking at the control propagator associated to H c (t) as the basic ob- 
ject for control design, 

U c {t) =Texpj-i^ dxH c (x)^ , (3) 

where T denotes time ordering and units such that h = 1 have been 
chosen. The analysis and complexity of the control problem are in- 
fluenced by various factors. Beside the available knowledge about the 
target system and the relevant noise generators S a , a critical role is 
played by the assumptions on the accessible control resources - in par- 
ticular, bounded vs unbounded control strengths and rates, and/or 
faulty control operations as opposed to perfect ones. 
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3 Bang-bang dynamical decoupling 



3.1 The simplest setting 

One can gain a concrete feeling on how the above strategy works by 
revisiting the simplest decoherence control scenario, namely a single 
spin-1/2 system (a qubit) diagonally coupled to an harmonic reservoir, 

H = Lu a z <g> l E + l s ® ^2 Vkblbk + <J Z ® ^ g k {b\ + b k ) , (4) 

fc fc 
where LOo,w k ,g k are real parameters, <j z = diag[+l, —1] is the z-Pauli 
matrix, and b\,b k are bosonic creation and annihilation operators. 
Since the interaction term Hse has opposite sign depending on whether 
the spin is up (|0)) or down (|1)), one may intuitively expect that it 
should be possible to average such contributions to zero by flipping the 
qubit faster than the typical time it takes to the environment oscillators 
to appreciably evolve. 

Formalizing such intuition has led to the first example of a bang- 
bang (BB) decoupling protocol |lf>j . In this case, the appropriate con- 
trol action is implemented by a train of resonant, infinitely short and 
strong (BB) 7r-pulses about the x (or y) axis, separated by a time inter- 
val At. In terms of the relevant control propagator U c (t), notice that 
this corresponds to a cyclic action with a period T c = 2Ai, that is, 
U C {T C ) = (n) 2 = Is, with 7T = i cxp(— ina x /2). The controlled evolu- 
tion of the qubit coherence element, poi(t), may be evaluated explicitly 
in this case. In the limit of arbitrarily fast control, this leads to the 
exact result 

lim p 01 (t = NT c )=p 01 (t = 0), (5) 

l c — >0,iv — >oo 

implying, in principle, the possibility to stroboscopically preserve quan- 
tum coherence indefinitely in time. In practice, the essential physical 
requirement making decoherence suppression possible is the ability to 
access control timescales T c over which the system-environment cou- 
pling can be coherently manipulated. Accordingly, if r c denotes the 
shortest correlation time associated with the reservoir dynamics, the 
condition T c <C t c is always necessary and sufficient in order to ap- 
proach the ideal limit of Eq. JSJ where decoherence is completely in- 
hibited. For the linear spin-boson model of Eq. Q, r c is essentially 
determined by the inverse of the highest frequency component present 
in the oscillator bath. 

By analogy with the well-know Carr-Purccll (CP) 7r-pulse sequence 
which is routinely used in NMR to remove unwanted phase evolu- 
tions , I will refer to this simple singlc-qubit decoupling protocol as 
CP-decoupling (see also for an explicit comparison between stan- 
dard and time-symmetric CP-schemes in this model and for the 
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extension to nonlinear spin-boson couplings). Interestingly, a proof- 
of-principle demonstration of CP-decoupling to suppress single-photon 
decoherence in an unbalanced Michelson polarization interferometer 
has been reported in j^J. Variants of this basic scheme (including 
27r-pulse protocols that are related to a z -decoupling) have found ap- 
plication in problems as different as suppression of collisional decoher- 
ence |32| . inhibition of spontaneous emission |33| . and error control in 
quantum search algorithms via atomic arrays |34j . 

3.2 The general case 

The above example can be generalized to BB decoupling of a generic 
(finite-dimensional) open quantum system S coupled to its environ- 
ment E as in Eq. JQ) ^|[5U|. Under the assumption that, as stip- 
ulated so far, Hq is time- independent and the applied control action 
is cyclic with a period T c , the stroboscopic controlled evolution of the 
combined system may be described by a propagator 

U(t = NT C ) =exp(-iH eff t) , NeN, (6) 

for a time-independent effective Hamiltonian H e f f . If, in addition, T c 
is sufficiently short, H e ff is accurately represented by the following 
lowest-order average Hamiltonian: 

H eff - H {0) = 1 f ° dx Ut (x)H U c (x) + 0(T C ) . (7) 

The term 0(T C ) accounts for higher-order corrections which can be sys- 
tematically evaluated through the appropriate Magnus expansion |28j , 
and converge to zero as the fast control limit T c — > is approached. 

The basic idea of a BB protocol is to appropriately map the time 
average contained in Eq. (JJJ into a group-theoretic average, by assign- 
ing the values of the control propagator U c (t) according to a discrete 
decoupling group Q, faithfully and unitarily represented on the state 
space Tis of S. Q is determined by the set of attainable BB control 
operations. If Q = {go, gi, ■ ■ ■ , <7l}j with L = \Q\ — 1, then T c is divided 
in L — 1 equal subintervals of length At, and BB decoupling according 
to Q is implemented by sequentially sampling the control propagator 

u c {t) from g, 

U c (t = (l-i)At + a)=g t -i, t=l,...,L, sG [0,At). (8) 

Thus, the control propagator is instantaneously changed from ge-i to 
ge through the application of a BB pulse pi = gtg\_ x at the end of 
each control subintcrval, tt = (-At. In this setting, CP decoupling 
corresponds to decoupling according to the simplest nontrivial group 
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J12, realized as Gcp = {I, o~ x } on TCs = C 2 , with BB control pulses 

Pi = P2 = <?x = 7T- 

The advantage of the group-theoretic framework is that it allows 
to straightforwardly construct controlled evolutions with well-defined 
symmetry properties with respect to Q, which can in turn be exploited 
to filter out the unwanted decohering terms. In particular, it was shown 
in [171 that the effective Hamiltonian of Eq. J7J) takes the form 

# (0) = w\ 2 g ? Ho9i ■ (9) 

11 ateg 

Because the above Hamiltonian is invariant under Q , [gi, H^] = for 
all £, this implies a controlled dynamics which is effectively symmetrized 
according to Q over any timcscalc longer than the averaging period T c . 
The exploitation of the symmetry structure enforced in this way by the 
controller is central to all aspects of the BB decoupling program, from 
the design of averaging schemes for relevant classes of environmental 
couplings |17| . to the characterization of universal procedures for de- 
coupled control j23], and the identification of dynamically generated 
noise-protected subsystems jSSJ- 

While being conceptually attractive thanks to the simplicity of the 
control design, BB decoupling schemes suffer from severe shortcomings 
in terms of their practical viability. In particular, the required control 
resources appear very unrealistic for two primary reasons: 

• Timescale requirements. Reaching control timescales T c as short 
as the minimum correlation time of the dynamics to be removed 
may imply prohibitively fast timescales for realistic open quan- 
tum systems. 

• Amplitude requirements. Even in the case where the control 
period T c is allowed to be reasonably long, instantaneous decou- 
pling pulses imply unbounded control strengths which are never 
affordable in practice. 

Thus, the relevant question is: Can we decouple under more real- 
istic assumptions? 

4 Relaxing amplitude requirements: Eu- 
lerian control design 

It turns out that the amplitude constraint can always be relaxed, at 
the expenses of making the control design slightly more sophisticated. 
The basic idea, introduced in [201, is to replace the piecewise constant 
BB propagator given in Eq. JHJ with a continuous control propagator 
able to ensure the same averaging, possibly with a reasonable overhead 
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Figure 2: Cay Icy graph of the group Z2 = {0,1}, represented as Q = {I,<r x } 
on C 2 . Points along the path which are specularly located relative to the center 
of the graph are correlated from the point of view of the averaging process. The 
parameter s <G [0, At) measures the elapsed time within each subinterval. 

in the required cycle time. This is accomplished by constructing the 
new control propagator U c (t) according to so-called Eulerian cycles 
on the Cayley graph G of Q that is, closed paths on G which possess 
the special property of using each edge of the graph once and once 
only [HfTj . 

The working principles of Eulerian control are best appreciated by 
comparing Eulerian and BB implementations in the simplest case of 
decoherence suppression in a single qubit already considered in Section 
3.1. In this case, Q = Qcp = {I, <?x}, and the general BB prescription 
of Eq. ijHJ takes the explicit form 



Physically, the first 7r-pulse may be thought of as instantaneously re- 
versing the sign of the underlying a z Hamiltonian, causing the evolu- 
tion during the second control subinterval to effectively retrace itself 
and lead to a zero average over T c . For Eulerian implementation, con- 
trol actions are no longer instantaneous, but rather distributed along 
the whole duration of each control subinterval At. The relevant Cay- 
ley graph for Q ~ Z 2 is depicted in figure 2. Because there is only one 
generator, 7 = a Xl the Eulerian prescription turns out to involve no 
time overhead with respect to the BB case, thus T c = 2At as before. 
Let h x (t) = f(t)cr x be any Hamiltonian realizing the generator 7 over 
At that is, 



u x(t) = Texp <^— i J dsh x (s) j , such that u x (At) = 7 = a x . (11) 



On the graph pictured in figure 2, an Eulerian cycle Ve is described 
by the sequence of edges Ve = (7,7)- Because the path determines 




- (a z + o x o z o x ) = - (a z - a z ) = . 



(10) 
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the sequence of control Hamiltonians to be turned on along the cycle, 
in this case one simply obtains a control propagator given by 



By inserting the above expression in Eq. (J7J), the net phase evolution 
over a cycle becomes 



By comparing Eq. I|13|) with Eq. (|10|) . one sees that the Eulerian 
protocol still achieves averaging by inducing an effective time reversal 
in the second control subinterval; however, such a reversal is now local 
in time, the evolution in the second interval continuously undoing the 
one in the first owing to the cancellations from opposite points in the 
graph (see figure 2). 

The generalization of the above procedure to generic decoupling 
protocols is described in [25] ■ Under mild assumptions on the relevent 
control Hamiltonians, one can always ensure that the same CJ-symmetri- 
zation of the BB limit is retained. However, Eulerian design turns out 
to be superior in two important respects: not only can any desired de- 
coupling protocol be fully implemented using bounded- strength Hamil- 
tonians, but the resulting schemes are also intrinsically stable against 
a large class of systematic control errors, thereby leading to enhanced 
robustness. 

5 Relaxing timescale requirements: Decou- 
pling of low-frequency noise 

Going back to reconsidering decoupling timescale requirements, it turns 
out that they can be substantially weakened for a wide class of envi- 
ronments generating low-frequency noise. I will focus here on the most 
notable representative instance, provided by noise processes exhibit- 
ing a 1// spectrum. Noise phenomena consistent with a 1// behavior 
have been identified to play a role in an astonishing variety of dynami- 
cal systems, ranging from physics to biology, economics, and more. At 
the quantum level, 1// noise due to a variety of material-specific fluc- 
tuation mechanisms is ubiquitous in solid-state devices, with a severe 




(12) 
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impact on the performance of both nanoelectronic circuits |37j and 
qubit realizations. In particular, 1/f noise due to fluctuating back- 
ground charges is currently regarded as one of the primary decoherence 
sources for superconducting charge qubits 

An appropriate schematization of 1/f noise effects requires going 
beyond the notion of a continuum of weakly coupled environmen- 
tal modes at equilibrium which is embodied by standard harmonic- 
oscillator reservoirs. A quantum-mechanical model capturing the dis- 
tinctive features of a realistic discrete environments generating 1/f 
noise was recently investigated in |39| . Within a semiclassical ap- 
proach that will suffice to the present discussion, a simple picture of 
1/f noise may be obtained by considering an incoherent superposition 
of relaxation processes from an ensemble of M independent bistable 
fluctuators. Let £,k{t) denote a random telegraph noise (RTN) signal 
describing the switching of an individual fluctuator between two values 
±Wfc/2 with a total switching rate 7^, and a corresponding Lorcntzian 
power spectrum Sfc(cj) = 7^/(7^ + cj 2 ), lo e E. Provided that the 
switching rates in the ensemble are assumed to be distributed in a 
range [y m i n ,lmax) with probability 
const 

p (7) = , const = [ln(7 roax /7 mi „)]" 1 , (14) 

7 

the total fluctuation S(t) = J^ fe (t) exhibits a power spectrum S(u>) = 
Y^k Sk{u) = A/\cj\, A > 0, in a frequency interval j min < w < "f ma x- 
For distributions of coupling strengths {vk} sufficiently peaked around 
their mean value (v) , the parameter A is directly related to the product 
Nd(v) 2 , rid = M/ logiptmax / Imin) being the number of fluctuators per 
noise decade. 

The simplest scenario for decoupling of 1/f noise is obtained by 
considering a single qubit in the presence of the semiclassical RTN 
disturbance. Thus, the target system is described by a time- dependent 
effective Hamiltonian of the form 

H (t) = H Q + H RTN (t) = [Sir, + Aa x ] + ~{t)a z , (15) 

for real parameters f2, A. While an appropriate extension of the basic 
average Hamiltonian formalism is needed in general to handle time- 
dependent target dynamics, the required control action can be imple- 
mented in this case through the same CP-protocol discussed for the 
spin-boson example of Section 3.1, either in the BB limit or in the Eu- 
lerian form illustrated in the previous section. A more detailed analysis 
of this problem is available in gOl (see also 02 SHI EH for related 
work). In what follows I only focus on some illustrative results in the 
pure dephasing limit, A = 0. In this regime, the stroboscopic evolution 
of the qubit coherence element may be exactly evaluated, 

|poi(* = NT C )\ = exp(-r c (i)) \p 01 (t = 0)| , (16) 
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time (1/fi) 



Figure 3: Suppression of 1// pure dephasing, fi = 1,A = 0, starting from the 
initial qubit state \ip) = (|0) + \l))/^/2. Solid lines: free evolution. For controlled 
evolutions, stroboscopic data (circles and diamonds) are connected by eye-guiding 
lines. In all cases, ~f ma x = 100. From top to bottom, the other spectral parameters 
are as follows: (a) -f min = 10~ 4 , (v) = 10~ 4 ; (b) j mm = 10~ 6 , (v) = 1(T 4 ; (c) 
Imin = 10~ 4 , (v) = 1Q~ 2 . Control parameters arc: (a) and (b) At = 1000 (dot 
dashed), At = 100 (dashed), At = 10 (dotted); (c) At = 10 (dot dashed), At = 1 
(dashed), At = 0.1 (dotted). 

the analytic expression of the controlled decoherence functional being 
given in 001- A formal result similar to Eq. (0) can be explicitly 
established in the continuous limit where T c — > 0, ./V — -> oo. Thus, com- 
plete compensation of 1// effects requires access to control timcscales 
At < 1/jmax, so as to average out the influence of the fastest fiuctua- 
tor present in the ensemble - as expected on general control-theoretic 
grounds. However, substantially slower decoupling rates might turn 
out to suffice in this case provided that the requirement of perfect 
suppression is lifted. 

The situation is pictorially summarized in figure 3. The three spec- 
tra differ due to the more or less pronounced Gaussian character of 
the underlying fluctuations, purely Gaussian dephasing (top plot) be- 
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ing closest to a formal representation in terms of a standard oscillator 
bath. As evidenced from the data, a 60% coherence recovery is always 
achieved by At values which are at least three orders of magnitude 
longer than expected from the inverse upper cutoff, the latter figure 
improving to five orders of magnitude for slower dephasing processes as 
in the top and middle panels. While such a favorable scaling in terms of 
decoupling resources might seem in a sense to naturally follow from the 
predominance of low-frequency modes in the noise profile, it is worth 
stressing that similar conclusions are not immediately applicable to 
the full class of possible 1// dynamics, the actual decoupling perfor- 
mance being in practice fairly sensitive to both the presence of strongly 
non-Gaussian noise sources and the operating point of the qubit |40| . 
Yet, this example convincingly shows the existence of realistic decoher- 
ence scenarios where decoupling might be achievable under affordable 
control resources. 

6 Conclusions 

Dynamical decoupling techniques offer a well-defined conceptual frame- 
work for addressing a variety of open-loop coherent-control problems 
for quantum systems. While their use as a tool for decoherence con- 
trol is particularly relevant, the applicability of decoupling methods 
and concepts extends to a broader setting, encompassing schemes for 
quantum-dynamical engineering on both physical and encoded degrees 
of freedom 

From a practical standpoint, the implementation of decoupling 
schemes poses stringent requirements on the needed control opera- 
tions. Although these are still beyond the reach of current experi- 
mental capabilities for many systems, progress in ultrafast coherent 
control is steady |45) . hopefully making decoupling implementations 
less demanding for such systems in a near future. On the theoretical 
front, both the notion of Eulerian control and the decoupling of 1// 
noise indicate how conclusions based on the simplest decoupling proto- 
cols or on general control-theoretic bounds may be overly pessimistic, 
and that significant improvement may still be possible in relevant situ- 
ations. It is my expectation that the convergence of better theoretical 
modeling and analysis, together with practical advances in quantum 
technologies, will improve the prospects that dynamical decoupling 
eventually becomes a method of choice for decoherence control. 
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